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I try to explain some recent progress in understanding the non-perturbative dynamics of hot non- 
Abelian gauge theories. The non-perturbative physics is due to soft spatial momenta |p| ~ g 2 T 
where g is the gauge coupling and T is the temperature. An effective theory for the soft field 
modes is obtained by integrating out the field modes with momenta of order T and of order gT in 
a leading logarithmic approximation. In this effective theory the time evolution of the soft fields 
is determined by a local Langevin-type equation. This effective theory determines the parametric 
form of the rate for hot electroweak baryon number violation as T = ng 10 log(l/(/)T 4 . The non- 
perturbative coefficient n is independent of the gauge coupling and it can be computed by solving 
the effective equations of motion on a lattice. 



I. INTRODUCTION 

Let me start by formulating the problem I am going 
to discuss: How can one calculate thermal expectation 
values like 



C(h - t 2 ) = (0(h)0(t 2 )) 



(1) 



in a non-Abelian gauge theory, when the leading order 
contribution is due to spatial momenta of order g 2 Tl The 
operator 0(t) is a gauge invariant function of the gauge 
fields Afj,(t,x.) at time t. When I said the leading order 
contribution is due to momenta of order g 2 T ', I referred 
to the Fourier components of the gauge fields entering 
0(t). Finally 



tr {(. ■ -)e- H ' T } 
ti{e~ H / T } 



(2) 



denotes the thermal average, where H is the Hamilto- 
nian. As one may anticipate, it is not possible to compute 
such a correlation function in perturbation theory. 

This problem arises in the context of electroweak 
baryon number violation at temperatures, well above 
lOOGeV, when the electroweak symmetry is unbroken. 
Then the SU(2) gauge theory is very similar to hot QCD. 
Fortunately, it is a bit simpler because the gauge coupling 
is small. 

In my talk I will try to explain how such correlation 
functions can be computed at leading order in the gauge 
coupling M. We will see that this requires the use of an 



effective theory which results from integrating out field 
modes with hard (p ~ T) and semi-hard (p ~ gT) spatial 
momenta [J Due to the smallness of the gauge coupling 
this can be done in perturbation theory. This effective 
theory turns out to have a relatively simple structure 
and it can be easily used for numerical calculations on a 
lattice. 

Most of the discussion will be restricted to a pure gauge 
theory, i.e., without matter fields. I will comment on the 
role of matter fields when necessary. 



II. HOT ELECTROWEAK BARYON NUMBER 
VIOLATION 

One of the most intriguing aspects of the electroweak 
theory is its non-trivial vacuum structure which is related 
to baryon number violating processes. These processes 
play an important role in understanding the observed 
baryon-asymmetry of the universe Q . 

Baryon number is not conserved in the electroweak 
theory due to the chiral anomaly. The baryon number 
current satisfies f\ 
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(3) 



where F is the SU(2) field strength tensor. The rhs of 
Eq. (0) is a total derivative, 



tr 



(4) 
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'For spatial vectors I use the notation k = |k|. Four-vectors 
arc denoted by if" = (fc°,k) and I use the metric K = 
k 2 -k 2 . 

''Baryon plus lepton number is violated in the electroweak 
theory while baryon minus lepton number is conserved. 
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FIG. 1. One dimensional slice of the gauge-Higgs field 
configuration space. The vertical axis is the energy of a 
field configuration which has the minimal energy for given 
Chern-Simons winding number. 

Thus, integrating Eq. (j^) over 3-space and over time from 
ti to tf we can write the change of baryon number B as 

B(t f ) - B(U) = n f (N CS (t f ) - NasiU)) (5) 

where 

Ncs = 3^ €ijk I d3x ( f z ai - r abCA i A " A t) (e) 

is the so called Chern-Simons winding number of the 
gauge fields. Unlike Nqs itself, the difference AiVcs = 
Ncs{tf) — Ncs(ti) is gauge invariant since it can be writ- 
ten as a space-time integral of FF. 

The Chern-Simons number is intimately related to the 
non-trivial topology of the gauge-Higgs field configura- 
tion space. There are paths in this space leading from the 
vacuum to the vacuum which cannot be continuously de- 
formed to a single point. The minimal field energy along 
such a path cannot be made arbitrarily small by con- 
tinuous deformation. There is a minimal energy barrier 
which has to be crossed. The corresponding field config- 
uration is the so called sphaleron. Following such a non- 
contractible path, the Chern-Simons number changes by 
an integer. 

The simplest mechanical analogue of this situation I 
can think of is a rigid pendulum in a gravitational field 
which can swing all the way around in a circle. The 
Chern-Simons number corresponds to 4>/(2ir) where <j) is 
the rotation angle. The vacuum corresponds to the pen- 
dulum at rest and a non-contractible loop corresponds to 
net rotation of the pendulum by an integer times 2ir. 

The gauge-Higgs field configuration space is of course 
much more complicated, it is infinite dimensional rather 
than one dimensional as in the case of the pendulum. 
Fig. [j] displays a one dimensional slice of this space. It 
corresponds to a path on which the field configuration 
for a given value of Ncs has a minimal energy. 

If there were only gauge and Higgs fields, nothing spec- 
tacular would happen in a vacuum to vacuum transi- 
tion which changes Ncs- However, in the presence of 



fermions, the change of Ncs is accompanied by a change 
of baryon number due to Eq. (^|). 

At zero temperature, iVcs-changing transitions are 
possible only by tunnelling. The tunnelling probability 
is so small that it has no measurable physical effect. 

This situation changes at high temperature. Then it 
becomes possible to make such transitions by thermal ac- 
tivation across the energy barrier. The suppression factor 
is then given by a Boltzmann factor rather than the tun- 
nelling amplitude. For T of the order of the electroweak 
phase transition or crossover temperature T c ~ lOOGeV, 
the transition probability becomes unsuppressed. 

When the temperature is sufficiently below T c , most 
transitions are passing the energy barrier close to the 
sphaleron and the transition rate T can be calculated in 
a saddle point approximation. 

For T ~ T c this approximation breaks down and for 
T > T c there is no sphaleron solution at all. Therefore 
one has to employ non-perturbative methods to calculate 
the rate for this case. It has become popular to call the 
rate for T > T c the hot sphaleron rate. The idea how to 
compute it is the following: The Chern-Simons number 
is expected to perform a random walk. Then, for large 
t, thermal average of [iVcs(i) — -^cs(0)] 2 grows linearly 
with time, 

([JVcs(t) - N CS (0)} 2 ) -» VtT (7) 

where V is the space volume. The coefficient T is the 
probability for a iVcs-changing transition per unit time 
and unit volume. The rate for baryon number violation 
Tab is proportional to T. Thus Tab can be obtained 
by computing the real time correlation function At 
very high temperature T>T C the thermal mass of the 
Higgs field becomes large so that it decouples. Then it is 
sufficient to evaluate Eq. ([?]) in the pure gauge theory. 

For the following discussion it is important to iden- 
tify both the relevant length scale R of the problem and 
the relevant size of the gauge field fluctuations A A ||. 
First of all, we have to consider fluctuations which are 
large enough to make a change ANcs of order unity. To 
estimate ANcs, recall that FF = EB. For the elec- 
tric fields we have E = — A in Aq = gauge. Thus 
J dt~E ~ AA. Note that the relevant time scale drops out 
in this estimate. Furthermore, we estimate B ~ AA/R 
and J d 3 x ~ R 3 . Thus we need 

ANcs ~ g 2 J d 4 x -E-B ~ g 2 R 2 (A A) 2 ~ 1. (8) 

In addition, we have to require that the energy of the 
relevant fluctuations does not exceed the temperature. 
Otherwise the transition rate would be Boltzmann sup- 
pressed. Proceeding as above we require 

energy ~ J d 3 xB 2 ~ R (AA) 2 < T. (9) 

Combining Eqs. (||) and (g) we find 
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R>(g 2 TY 



AA ~ (gR)- 



(10) 



On the other hand, we know that the correlation length 
of magnetic fields in a hot non-Abelian plasma is of or- 
der ig 2 T)~ x . Thus we can have only J?<(g 2 T) _1 which 
finally gives 



R ~ (g 2 TY 



AA ~ gT. 



(11) 



As I stated above, the relevant time scale t, which 
we need to know in order to estimate the hot sphaleron 
rate, has dropped out in these considerations. It has to 
be determined from the dynamics and it cannot be ob- 
tained from purely thermodynamic considerations. Once 
we know t we can estimate the hot sphaleron rate as 



r - t^R- 3 . 



(12) 



While Eq. ( |ll| ) has been well established for a quite a 
while, the characteristic time scale has been understood 
only recently jy. I hope that, by the end of my talk, you 
will be convinced that the correct answer is given by Eq. 



III. NON-PERTURBATIVE PHYSICS AT HIGH 
TEMPERATURE 

After the discussion of electroweak baryon number vio- 
lation we now return to our original problem, which is to 
calculate correlation functions like ([I]). 

In the previous section we have seen that the rele- 
vant length scale for electroweak baryon number viola- 
tion is (<? 2 T) _1 corresponding to momenta p ~ g 2 T. It is 
well known that finite temperature perturbation theory 
in non-Abelian gauge theories breaks down for momenta 
as small as g 2 T One can see this by considering 

the typical size of the transverse gauge field fluctuations 
which can be estimated from the equal time correlation 
function. For the transverse gauge fields we have 



A t (x) ~ gT. 



(13) 



Thus the two terms in the covariant derivative di — igAi 
are of the same size which makes perturbation theory 
impossible. Note that the thermal fluctuations are of the 
same size as the change AA which we have estimated in 
Eq.©. 

When I said that perturbation theory does not work 
for soft momenta, one may ask: What about the plasmon 
at rest which has p = and the celebrated calculation 
of it's damping rate H? The point is that the plasmon 
oscillation has an amplitude ~ g 2 T , i.e., it is a small os- 
cillation. What we are interested in is the time evolution 
which makes the field change by an amount equal to it's 
typical size ~ gT. 

For the case of equal time correlation functions the 
non-perturbative physics associated with the scale g 2 T 



is relatively well understood. It is determined by an ef- 
fective 3-dimensional pure gauge theory. This effective 
theory is the result of integrating out the field modes 
with non-zero Matsubara frequency ( "dimensional reduc- 
tion" ) and the 0-component of the gauge fields . For 
sufficiently small coupling this can be done in perturba- 
tion theory. The 3-dimensional theory can be used for 
3-dimensional Euclidean lattice simulations which allow 
for a much higher precision than 4-dimensional ones ||. 
It should be noted that the use of dimensional reduction 
is more a matter of practical convenience rather than a 
matter of principle. The full 4-dimcnsional can also be 
evaluated directly in a 4-dimensional Euclidean lattice 
simulation. 

The situation is different for unequal time correlation 
functions. Dimensional reduction is not possible be- 
cause we have to consider time dependent quantities. 4- 
dimensional Euclidean lattice simulations are of no use 
either: Formally, real time correlation functions can be 
obtained by an analytic continuation of their imaginary 
time counterparts. However, it is not possible to do an 
analytic continuation of a function which can be evalu- 
ated only numerically. 

One simplification of the problem is that the soft field 
modes behave classically. This is due to the fact that 
the momentum scale of interest is small compared to the 
temperature. Then the number of field quanta in one 
mode with wave vector p, given by the Bose-Einstein 
distribution function 



nip) 
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e p/T _ i 



T 
P 



(14) 



is large. In this case we are close to the classical field 
limit. Thus the dynamics of the soft field modes should 
be governed by classical equations of motion. |**| 

The first estimate of the hot sphaleron rate, due to 
Khlebnikov and Shaposhnikov ||, was based on the as- 
sumption that only the soft fields play a role while the 
high momentum modes decouple at leading order. Then 
the only scale in the problem would be g 2 T and, on di- 
mensional grounds, the hot sphaleron rate would have to 
be of the form 



r KS - K(g 2 T) 4 



(15) 



with a non-perturbative numerical coefficient n. Subse- 
quently attempts were made to compute k by solving the 
classical Hamiltonian equations of motion for the gauge 
fields and measure T. 

However, it turns out that in the case of unequal time 
correlation functions is much more complex than in the 
equal time case. The main difficulty is to understand the 
role of the different momentum scales which are present 



*Tt should be noted that this argument is somewhat hand- 
waving and it is not easy to justify quantitatively |l0|-|l2|]. 
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in the problem. We will see that both momenta of order 
T and momenta of order gT are relevant to the leading 
order dynamics of the soft fields. 

Therefore one has to integrate out the high momentum 
modes to obtain an effective classical theory. Schemati- 
cally, the equations of motion will be of the form 



SS oB [A] 
5A u (x) 







where the effective action S e s [A] consists of the tree level 
action and the terms which are generated by integrating 
out the field modes with p ~ T and p ~ gT. Deriving 
this equation of motion will be the main subject of my 
talk. 

The reason why the hard modes have to be integrated 
out to obtain an effective classical theory for the soft 
modes is pretty obvious: For p ~ T the Bose-Einstein 
distribution is of order 1 so that these modes certainly 
do not behave classically. For semi-hard momenta, how- 
ever, we have n{p) <C 1. The main reason why they have 
to be integrated out as well is that it does not seem to be 
possible to construct a lattice formulation of the effective 
theory for momenta of order gT and g 2 T which allows to 
take the continuum limit (cf. Ref. p3|). Another reason 
is that after integrating out the scale gT we will obtain 
an effective theory for the soft field modes which con- 
tains only one length scale and one time scale. Then it 
will be trivial to estimate the parametric form of the hot 
sphaleron rate. 



IV. INTEGRATING OUT THE HARD MODES 

The hard modes constitute the bulk of degrees of free- 
dom in the hot plasma. Their physics is that of almost 
free massless particles moving on straight lines. Nev- 
ertheless, they have a significant influence on the soft 
dynamics because they are so numerous. 

Integrating out the hard modes means that we have to 
calculate loop diagrams with external momenta <C T and 
internal momenta of order T. This generates effective 
propagators and vertices for the field modes with p <C 
T. To leading order, we can restrict ourselves to one 
loop diagrams and we can make a high energy or eikonal 
approximation of the integrand. The result is nothing 
but the so called hard thermal loops Q . The generating 
functional of the hard thermal loops has the remarkable 
property to be gauge invariant. 

The main difference between QED and non-Abelian 
theories is that for the former the only hard thermal loop 
is the 2-point function, while for the latter there are also 
hard thermal loop n-point functions for all n. As we will 
see later, this has a significant effect on the soft dynamics, 
it is in fact qualitatively different in Abelian and non- 
Abelian theories. 

For the moment, however, we restrict the discussion to 
the hard thermal loop 2-point function. This will make 



clear that the hard modes have a significant influence on 
the soft dynamics which was pointed out by Arnold, Son, 
and Yaffe @. 

We have to consider the transverse (or magnetic) com- 
ponents of the gauge fields. The transverse polarisation 
function is given by 



5Ht(P) = X -m 



+ P0 1 



Pi 



dn v i 

4-7T V-P 



(16) 



where to|, is the leading order Debye mass squared. In a 
pure SU(iV) gauge theory, it is given by 



-Ng 2 T 2 . 
3 y 



(17) 



In the hot electroweak theory it receives additional con- 
tributions due to the Higgs and fermion fields. Further- 
more, = (l,v), and the integral J dfl v is over the 
directions of the unit vector v, |v| = 1 . 

For the naive "natural" scale of the problem, i.e., po, 
p ~ g 2 T we have 



6Tl t (P)~g 2 T 2 ( Po ~g 2 T, p~g 2 T), 



(18) 



i.e., (5II t is much bigger than the tree level kinetic term. 
Therefore 8H t must be resummed. The resulting trans- 
verse propagator 



*A t (P) 



1 



-P 2 + sn t (p) 



(19) 



is smaller than the tree level propagator by a factor g 2 
when po ~ S 2 ^P P ~ 9 2 T '• That means that for the 
naive "natural" frequency there are only small perturba- 
tive fluctuations which are of the same amplitude as the 
plasmon oscillation. It also means that the estimate ( |l5| ) 
cannot be correct. 

In order to obtain large non-perturbative fluctuations 
we obviously need to make <5n t smaller. We know that 
Silt vanishes PI for po — > 0. Thus we have to consider 
small values of po . In the limit po <C p we have 



fllt(P) 



2 PO 
P 



(20) 



which is purely imaginary. The imaginary part of <5n t (P) 
reflects the fact that the dynamics of the soft modes is 
Landau-damped by the hard modes. In fact, it is over- 
damping because the damping term is much larger than 
the tree level kinetic term p 2 . In the same limit the trans- 
verse propagator (nm becomes 



*At(P) 



1 



p^ po + l^p 



(21) 



^The longitudinal hard thermal loop polarisation tensor 
stays large even for po — * 0. Thus only the transverse fields 
can produce large, non-perturbative fluctuations. 
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with 



7 P 



4p3 



(22) 



Now we see how small po must be in order to obtain the 
desired large fluctuations, namely 



corresponding to the time scale 



(23) 



(24) 



Based on these considerations Arnold, Son and Yaffe es- 
timated the hot sphaleron rate as 



r ASY = ^ 10 r 4 



(25) 



Subsequently, attempts were made by Huet, Son [ pXi|Jl7[ 
and Arnold Jl8| to construct a numerical algorithm to 
compute the coefficient k on a lattice. 

In this discussion it is assumed that interactions do 
not change these order of magnitude estimates. As we 
will see in the next section, this is not the case: The field 
modes with momenta of order gT lead to interactions of 
the hard modes which cannot be neglected relative to the 
hard thermal loops. 



V. BEYOND HARD THERMAL LOOPS 

I will now argue that the hard thermal loop approxi- 
mation for the high momentum modes (p 3> g 2 T) is not 
sufficient to obtain the correct effective theory for the 
soft modes. 

Consider a hard thermal loop and imagine adding a 
self energy insertion on an internal line. This gives the 
diagram in Fig. [|a in which the hard loop momentum Q is 
on shell, Q 2 = 0. Since the external momentum P is soft, 
the momentum Q + P is almost on shell. The order of 
magnitude of the self energy insertion is well known from 
the calculation of the lifetime of hard particles [[l4| . From 
loop momenta K of order gT one gets a contribution 
~ g 2 T 2 . Compared to the hard thermal loop, diagram 
^a also contains an additional propagator which can be 
approximated as 



1 



1 1 



1 1 



(Q + P) 2 2qv-P T v-P' 
Thus we can estimate the diagram in Fig. Ha as 



nto(p) ~(Sn t (p) 



v-P 



(26) 



(27) 



and, by power counting, this diagram is not suppressed 
relative to the hard thermal loop ( |l6| ) when P ~ g 2 T. 
Note that in this estimate it did not play any role that 
we are considering a non-Abelian theory. The same type 




Q 



K 



(a) 



0>) 



FIG. 2. Two loop contribution to the polarisation tensor 
for soft external momentum P (po, p < g 2 T). The momentum 
Q is hard (go, q ~ T). The thick lines denote propagators 
with momenta of order T. The thin lines are hard thermal 
loop resummed gauge field propagators carrying momentum 
K with fc , k ~ gT. 



of diagram is also present in QED and the estimate would 
be exactly the same. However, there is an essential dif- 
ference between the Abelian and the non-Abelian case 
which will become clear in a moment. 

Similarly, one can see that the ladder-type diagram ||d 
is as large as the hard thermal loop as well. Again, this 
estimates are the same in QED and QCD. However, in 
QED the large contributions from the two diagrams turn 
out to be the same with opposite signs and cancel. 

I will now argue that this cancellation does not oc- 
cur in a non-Abelian theory. Let me first give a formal 
exp lanation in terms of diagrams. At the end of Sect. 
VII, I will come back to this point and I will give a more 
intuitive argument. 

Here is the diagrammatic argument: We have seen that 
the two diagrams both contain a hard and a semi-hard 
loop momentum. Imagine now that we calculate them 
in two steps: First do the integral over the hard mo- 
mentum with K kept fixed. This can be interpreted as 
follows: The first step generates an effective 4-point ver- 
tex for momenta of order gT and g 2 T. The result of such 
a calculation is well known: At leading order this is noth- 
ing but the hard thermal loop four point function. And 
we also know that there is no such vertex in QED, even 
though individual diagrams would give such a contribu- 
tion. Summing over the permutations of external lines 
the large contributions cancel. In a non-Abelian theory 
each of these diagrams comes with a different colour fac- 
tor and the sum over permutations gives a non- vanishing 
result. In the second step we integrate over momenta of 
order gT which can be done using the well known expres- 
sion for the hard thermal loop 4-point function. 



VI. INTEGRATING OUT THE SEMI-HARD 
MODES 

We have seen in the previous section that, in order to 
obtain an effective theory for the soft momentum fields, 
we cannot restrict ourselves to the hard thermal loops. 
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(a) 



(b) 



FIG. 3. One loop contributions to the polarisation tensor 
in the "hard thermal loop effective theory". The heavy dots 
are hard thermal loop vertices. Otherwise the notation is the 
same as in Fig. El 



Not only the hard modes but also the semi-hard ones af- 
fect the leading order dynamics of the soft fields. At first 
sight this appears like an additional complication. Some- 
what surprisingly, things become much simpler instead. 

We have also seen that it is convenient to proceed in 
two steps: In the first step we integrate out momenta 
of order T which yields the hard thermal loop effective 
theory. In the second step this effective theory is used to 
integrate out the scale gT. In this way the calculation 
will be simplified because we have a partial cancellation, 
corresponding to the QED-type contributions, already 
"built in" . 



compensated by one additional factor v-P ~ g 2 T in the 
denominator. Eq. (^7j) is even larger than ([If]) by a factor 
oflog( 5 T/ M ) 

Note that ( p9|) is not transverse, which would be nec- 
essary for the effective theory, which we are going to de- 
rive, to be gauge invariant. There is, however, another 
one loop diagram in the hard thermal l oop effective the- 
ory which is of the same size as Eq. ( |27|) and which is 
depicted in Fig. |^b. With the same approximations as 
above one obtains |f9|] 



x po 



dfi vi v-in fdQ V2 v 2v (vi-v 2 ) 2 



4tt v x -P J An v 2 -P ^\ _ ( Vr v 2 ) 



(30) 



One easily verifies that the sum of Eqs. ( |29| ) and ( |30| ) is 

transverse, 



p» (n$(p) + ii<$(p) 



= o. 



(31) 



B. Higher loops 



A. One loop diagrams 

We start by considering the one loop diagram in the 
hard thermal thermal loop effective theory depicted in 
Fig. ||a. It corresponds to the sum of the diagrams of the 
original theory in Fig. ||. We neglect the external soft 
momentum relative to the semi-hard one whenever it's 
possible. In this way we obtain an expression in which 
the loop integration is logarithmically infrared divergent. 
This should not bother us since we do not want to com- 
pute the diagram exactly. All we want to do is to in- 
tegrate out momenta of order gT. Therefore we need 
a scale /i which separates semi-hard momenta from soft 
ones, 



g 2 T « fj, « gT. 



(28) 



This scale appears as a lower limit in our integral so that 
it is no longer divergent. The integral is difficult evaluate. 
What it not so difficult, is to compute the terms which 
are singular for fi — > 0. One finds |l9| 



U^(P) = -^m 2 D Ng 2 T\og^ 



(29) 



This result is gauge fixing independent in a general co- 
variant gauge. 

Obviously, Eq. ( gjj) is not suppressed relative to the 
hard thermal loop selfenergy (Jig): The prefactor g 2 T is 



There are higher loop diagrams in the hard thermal 
loop effective theory which are not suppressed relative 
hard thermal loops either. Obviously these terms need 
to be resummed. As I already said, the use of the hard 
thermal loop effective theory has simplified the calcula- 
tion of ( p9| ) and ( |30| ) quite a bit. Nevertheless, higher 
loop diagrams are still difficult to calculate because the 
number of terms in the hard thermal loop rt-point func- 
tions grows rapidly with increasing n. 

It is much more efficient to use another formulation of 
the hard thermal loop effective theory, which is in terms 
of kinetic equations. The degrees of freedom in this de- 
scription are the gauge fields A fi (x) and, in addition, the 
fields W(x, v) which transform under the adjoint repre- 
sentation of the gauge group. The fields W(x, v) depend 
both on the space-time coordinates and on the unit vec- 
tor v. They describe the deviation of the distribution 
of hard particles with 3-velocity v from thermal equilib- 
rium. 

The equations of motion for these fields are the non- 
Abelian generalisation of the Vlasov equations for a QED 
plasma [E0|, 



[D fl ,F»»(x)] = 



47T 



■v"W(x,v), 



[v ■ D,W(x,v)] 



E(a:), 



(32) 



(33) 



where E is the non-Abelian electric field, and the 4- vector 
v is defined as in Eq. @. The rhs of Eq. (§|) is the 
current due to the hard particles. The conserved Hamil- 
tonian corresponding to Eqs. (|32|) and (|3|) is 2~jj22| 



G 



H = / d 3 xtr{ E(x) • E(x) + B(x) • B(x) 



4tt 



—v w(x,v), 



(39) 



4tt 



-W(x,v)W(x,v) 



(34) 



One advantage of this formulations that this theory is 
local. What is even more important, however, is that the 
algebraic complexity is significantly reduced. 

What does it mean to integrate out the semi-hard field 
modes in this formulation? First we have to split the 
fields into two pieces. The fields A, E and W are decom- 
posed into soft and a semi-hard modes Q 



A -> A + a 
E -> E + e 
W ^W + w. 



(35) 



The soft modes A, E and W contain the spatial Fourier 
components with p < fi while the semi-hard modes a, e 
and w consist of those with k > fj,. 

The interpretation of W and w is the following: Both 
describe deviation of the distribution of the hard particles 
from thermal equilibrium. W (w) is the slowly (rapidly) 
varying piece of this distribution varying on length scale 
greater (less) than 

After this split we have two sets of equations of motion, 
one for the soft fields and one for the semi-hard fields. 
Due to the non-linear terms these sets are coupled. 

The equations for the soft fields can be written as 



[D ll) F^(x)}=m 2 D J 



dn^ 

47T 



v v W{x,v), 



[v ■ D, W(x, v)] = v E(x) + £(x, v ), 



where 



Z a (x,v) = -gf abc (v-a b (x)w c (x,v)) 



soft 



(36) 
(37) 

(38) 



The subscript "soft" indicates that only spatial Fourier 
components with p < ji are included. The field strength 
tensor and the covariant derivatives in Eqs. (pfj) and ( |37| ) 
contain only the soft gauge fields. Note that terms which 
couple soft and semi-hard fields have been neglected in 
Eq.@. 

Since the semi-hard modes can be treated perturba- 
tively, one can neglect their self-interaction. One can 
also neglect interactions in the first Vlasov equation, so 
that we have 



**For notational simplicity I do not introduce new symbols 
for the soft modes. From now on A, E and W will always 
refer to the soft fields only. 

§ § The structure constants f abc are defined such that 



v-dw(x, v) — v-e(x) + h(x, v), 



(40) 



where / M!y = d^a? — d v . The only interaction which af- 
fects the time evolution of the semi-hard fields at leading 
order appears in Eq. (E(j), it is due to the term 



h a {x,v) 



-gf 



(the 



v-A b (x)w c (x,v) 



■va b {x)W c {x 1 v) 



(41) 



Integrating out the scale gT means that we solve the 
classical equations of motion ( |39| ) and ([lo]) for the semi- 
hard fields on the soft background. Then we have to 
perform the thermal average over their initial conditions 
using the Hamiltonian ( |34]) . In this way we eliminate the 
semi-hard fields and we end up with equations of motion 
for the soft fields only. This can be done in perturba- 
tion theory, where one expands in powers of h(x,v). We 
obtain an expansion 



a(x) = a (x) + ai(x) + a 2 (x) + ■■■ 
w(x, v) = w (x, v) + wi(x, v) + w 2 {x, v) + 



(42) 



in which the n-th order term contains n powers of the 
soft fields A, W. The terms do, wq are independent of 
the soft background and they only depend on the initial 
conditions for the semi-hard fields. Note that the solution 
( pl2| ) is still formal since it contains the yet unknown non- 
pcrturbative soft fields. 

Inserting ((42[) into Eq. (B8T) we obtain a series 



(43) 



Each term in ( }43| ) is bilinear in free fields ao and wq. 
The term £„ is of n-th order in the fields A and W. This 
expression has to be plugged into the second Vlasov equa- 
tion IB7l) for the soft fields, which then becomes 



(44) 



[vD,W(x,v)\ =v-E(x) 

+ €q(x,v) +£i(x,v) +£, 2 (x,v) + 



As it stands, Eq. (|44j) is not very useful yet. On the 
rhs we have an infinite series of terms, each depending on 
the yet unknown soft fields and on the initial conditions 
for the semi-hard fields. In general, the average over 
these initial conditions can be performed only after the 
equations for the soft fields have been solved. 

Fortunately, we can simplify Eq. ( [f4"| ) significantly since 
we are only interested in the leading order dynamics of 
the soft field modes. Recall that each term is bilinear 
in the fields ao and wq. Furthermore, ao and wq are linear 
in the initial values for their time evolution. After solving 
the equations of motion for the soft fields (^) and (H), 
one has to average over initial conditions. Due to the 
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scale separation the average of the bilinear terms of semi- 
hard fields can be approximated by disconnected parts. 
But this just means that we can perform the average over 
initial conditions for the semi-hard fields already in the 
equation of motion ([li])! The only exception is the term 
£o because the thermal average of this term vanishes due 
to colour conservation: Thermal averages like (d^dg) are 
proportional to S ab . Inserted into Eq. ( p8|) this gives zero 
due to the antisymmetry of the structure constants f abc . 
The highest correlation function of £o we have to take into 
account is the two point function (£o(xx,V\)£o(x2,V2)). 
Put differently, the term £0 acts like a Gaussian random 
force. It is random because it is independent of the soft 
fields. 

Since the expectation value of £0 vanishes, it is not suf- 
ficient to keep only this term on the rhs of Eq. . We 
also have to take into account £1 which we can approxi- 
mate by 



(45) 



where ((■•■)} denotes the average over initial conditions 
for the semi-hard fields. The higher order terms £„ in 
Eq. ( |44| ) can be neglected. Then we have 



[v-D,W(x,v)] ~ v-E(x) 



(46) 



in which the correlation functions of £0 can be treated as 
Gaussian. 



VII. THE LOGARITHMIC APPROXIMATION 

The next step towards the effective equations of mo- 
tion for the soft fields is to evaluate the terms in the 
second line of Eq. ([l6]) . As in the diagram calculation in 
Sec. VIA one encounters terms which are logarithmically 
sensitive to the separation scale /1. Only these logarith- 
mic terms will be kept. The reason why this is sufficient 
will become clear in a moment. 

We have seen that all we need to know about the noise 
term £0 is its 2-point function. The leading logarithmic 
result reads 



x5 ab S ( - 4 \x 1 -x 2 ), 



(47) 



with 



I(v,vi) = -5 {s2) (v- vi) + 



1 (vvr) 2 
tt 2 Vl-(v-vr) 2 



(48) 



where ) is the delta function on the two dimensional 
unit sphere: 



dO Vl /(vi)^ s2 )(A 



vi 



/(v) 



(49) 



With the same approximations the result for the last term 
in Eq. © is 



Ux(x,v))) = Ng 2 T log 



9l 



4?r 



-l{v 1 vx)W{x 1 v x ). 



(50) 



Inserting this into Eq. (46) we obtain 



[v-D,W(x,v)} =vE(i)+Jo(i,«i) 

"dfi 

M 



4tt 



/(l),t)l)^(x,Di). (51) 



I will now argue that, for the leading order dynamics of 
the soft fields, the lhs of Eq. @ can be neglected. The 
only spatial momentum scales which are left in the prob- 
lem are /i and g 2 T . The field modes we are ultimately 
interested in, are the ones which have only momenta of 
order g 2 T . The cutoff dependence on the rhs must drop 
out after solving the equations of motion for the fields 
with spatial momenta smaller than fx. Thus, after the /it- 
dependence has cancelled, the logarithm must turn into 
log(gT/(g 2 T)) = log(l/ 5 ). 

We will now simplify Eq. ( p3l| ) by neglecting terms 
which are suppressed by inverse powers of log(l/g). We 
introduce moments of W(x, v) and ^o(x, v), 



' 47T 



■v ln W(x,v), 
■v ln €o(x,v). 



(52) 



Taking the first moment of Eq. (^lj) gives 

1 



[D ,W*{x)]-[D>,W* 3 {x)\ 
Ng 2 T 



-E\x) 



in 



\og(l/ g)W l (x). 



(53) 



Assuming that (cf. Ref. fl9|l ) is of the same size as 
W l (x) and D <g 2 T, the lhs of Eq. Q is logarithmically 
suppressed relative to the term oc W l {x) on the rhs and 
can be neglected. Then we can trivially solve ( |53| ) for 

W l {x), 



W l {x) 



47T 



iV ff 2 Tlog(l/<?) \3 



(54) 



Inserting the current j i (x) — m^ i W' l {x) into the rhs of 
Eq. (H) gives 



[D,,F^(x)}^ 1 E\x)+C(x) 
where we have introduced 



47rmp 



Ng^T log(l/g) 



(55) 



(56) 
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and 



47rmQ 



Therefore the hot sphaleron rate, at leading order, has 
the form 



37V. 9 2Tlog(l/ ff )- 



(57) 



The term £ is a Gaussian white noise. Its correlator is 
easily obtained from Eq. (f47|), 

(C a (xiK Jb (x 2 )) = 2T 1 ^5 ab 5 {i \x l - x 2 ). (58) 

Eq. ( |55| ) is the main result of our calculation. It is re- 
markable that, at the order we are considering, the effect 
of the high momentum modes can simply be described 
by a local damping term, together with a Gaussian white 
noise which keeps the soft modes in thermal equilibrium. 
This surprising result is somewhat counter-intuitive since 
Landau damping is generally known as a non-local effect. 
Qualitatively, it can be understood as follows: Recall 
that the hard particles have a mean free path of order 
(<7 2 Tlog(l/<7)) _1 [Q. To understand the physical pic- 
ture behind the logarithmic approximation, we have to 
imagine that log(l/g) is a large number. Then the mean 
free path of the hard particles is small compared to the 
length scale (g 2 T) _1 on which the soft non-perturbative 
dynamics occurs. That is, the soft field modes are Lan- 
dau damped only over a small length scale. 

There is an essential difference compared to Abelian 
gauge theories. There the mean free path is of order 
(g 2 Tlog(l/g)) -1 as well. However, the mean free path is 
determined by the total cross section of the hard parti- 
cles which is dominated by small angle scattering. Thus, 
in a typical scattering event, the momentum of a charged 
particle in a QED plasma is hardly affected. The Landau 
damping can continue as a coherent process after such a 
scattering occurs. In a non-Abelian theory even a scat- 
tering under arbitrarily small angle can change the colour 
charge of a hard particle. It is precisely this colour ex- 
change which leads to a loss of coherence in the process 
of Landau damping. 



VIII. THE HOT SPHALERON RATE 

With Eq. (|59) we are now able to estimate the charac- 
teristic time scale t of non-perturbative gauge field fluc- 
tuations and thus of the rate for hot electroweak baryon 
number violation. Neglecting the term [Dq, F 0l (x)] (see 
below) the lhs can be estimated as 



[D 3 ,F Jl {x)\ ~ R- 2 AA ~ 5 4 T 2 AA 
On the rhs we have in Aq = gauge 

T AA 

7E ■ 



(59) 



(60) 



log(l/ff) * 
Putting @ and @ together we find 

t-O^loga/ffjT)- 1 . (61) 



r = Kff 10 log(l/ 5 )T 4 



(62) 



where k is a non-perturbative coefficient which does not 
depend on the gauge coupling. 



IX. AN EFFECTIVE THEORY FOR THE SOFT 
MODES 

I will now discuss how the non-perturbative coefficient 
k in Eq. ([62]), and more generally, correlation functions 
like (|l|) can be calculated on a lattice for which it is con- 
venient to work in the Aq = gauge. 

The time scale for non-perturbative dynamics is much 
larger than the corresponding length scale. Therefore 
time derivatives on the lhs of Eq. (pq) are negligible and 
the non-perturbative dynamics of the soft gauge fields at 
leading order is correctly described by the equation 



[D j ,F* i {x)]=->yA i (x) 



(63) 



which should be easy to implement in a lattice calcula- 
tion. 

After solving Eq. (|6^), the result has to be plugged 
into the operator 0(t) of interest. The corresponding 
correlation function is then given by the average over the 
random force £ using Eq. (pq). 



X. SUMMARY AND DISCUSSION 

We have obtained an effective theory for the non- 
perturbative dynamics of the soft field modes by integrat- 
ing out the hard (p ~ T) and semi-hard modes (p ~ gT) 
in perturbation theory. This effective theory is described 
by the Langevin equation ([33]). 

Furthermore, we have determined the parametric form 
of the hot electroweak baryon number violation rate at 
leading order. It contains a non-perturbative numerical 
coefficient which can be evaluated using Eq. (63). 

One would expect corrections to Eq. (p2) to be sup- 
pressed by a factor (logfl/g)) -1 , which, in the elec- 



troweak theory, is not small ***. It would therefore be 
interesting, both from the theoretical and from the prac- 
tical point of view, to see how sub-leading terms can be 
computed. 



*We have g ~ 0.66 for T ~ lOOGeV which gives log(l/<?) 



0.4. 
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Note added: While this paper was being completed, 
I received a preprint by Arnold, Son and Yaffe p3| , in 
which they derive Eq. (^3|) using the concept of colour 
conductivity. They also show that the effective theory 
described by this equation is insensitive to the ultravio- 
let which means that its lattice implementation does not 
require any renormalization. 
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